### Placebo tests for the two events with significant effects

## Original Lampedusa
# 90 days before
rdlamp1 <- RDestimate(factor ~ lamp, data = data61, cutpoint = -90, bw = 30)
summary(rdlamp1) 

rdlamp2 <- RDestimate(factor ~ lamp, data = data62, cutpoint = -90, bw = 30)
summary(rdlamp2)

# 75 days before
rdlamp1 <- RDestimate(factor ~ lamp, data = data61, cutpoint = -75, bw = 30)
summary(rdlamp1) 

rdlamp2 <- RDestimate(factor ~ lamp, data = data62, cutpoint = -75, bw = 30)
summary(rdlamp2)

# 60 days before
rdlamp1 <- RDestimate(factor ~ lamp, data = data61, cutpoint = -60, bw = 30)
summary(rdlamp1) 

rdlamp2 <- RDestimate(factor ~ lamp, data = data62, cutpoint = -60, bw = 30)
summary(rdlamp2)

# 45 days before
rdlamp1 <- RDestimate(factor ~ lamp, data = data61, cutpoint = -45, bw = 30)
summary(rdlamp1) 

rdlamp2 <- RDestimate(factor ~ lamp, data = data62, cutpoint = -45, bw = 30)
summary(rdlamp2)

# 30 days before
rdlamp1 <- RDestimate(factor ~ lamp, data = data61, cutpoint = -30, bw = 30)
summary(rdlamp1) 

rdlamp2 <- RDestimate(factor ~ lamp, data = data62, cutpoint = -30, bw = 30)
summary(rdlamp2)

# 365 days before
rdlamp1 <- RDestimate(factor ~ lamp, data = data61, cutpoint = -365, bw = 30)
summary(rdlamp1) 

rdlamp2 <- RDestimate(factor ~ lamp, data = data62, cutpoint = -365, bw = 30)
summary(rdlamp2)




## 23.03.2017
# 90 days before
rd10 <- RDestimate(factor ~ ten | cntry, data = data81, cutpoint = -90, bw = 30)
summary(rd10)

rd10b <- RDestimate(factor ~ ten | cntry, data = data82, cutpoint = -90, bw = 30)
summary(rd10b)

# 75 days before
rd10 <- RDestimate(factor ~ ten | cntry, data = data81, cutpoint = -75, bw = 30)
summary(rd10)

rd10b <- RDestimate(factor ~ ten | cntry, data = data82, cutpoint = -75, bw = 30)
summary(rd10b)

# 60 days before
rd10 <- RDestimate(factor ~ ten | cntry, data = data81, cutpoint = -60, bw = 30)
summary(rd10)

rd10b <- RDestimate(factor ~ ten | cntry, data = data82, cutpoint = -60, bw = 30)
summary(rd10b)

# 45 days before
rd10 <- RDestimate(factor ~ ten | cntry, data = data81, cutpoint = -45, bw = 30)
summary(rd10)

rd10b <- RDestimate(factor ~ ten | cntry, data = data82, cutpoint = -45, bw = 30)
summary(rd10b)

# 30 days before
rd10 <- RDestimate(factor ~ ten | cntry, data = data81, cutpoint = -30, bw = 30)
summary(rd10)

rd10b <- RDestimate(factor ~ ten | cntry, data = data82, cutpoint = -30, bw = 30)
summary(rd10b)

# 365 days before
## No data available



